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Abstract 

The pole mass of the gluino is diagrammatically calculated to the two-loop order 
in SUSY QCD as a function of the running parameters in the lagrangian, for the 
gluino and squarks sufficiently heavier than the quarks. The 0(aj?) correction shifts 
the gluino mass by typically 1-2 %, which may be larger than the expected accuracy 
of the mass determination at future colliders. 
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If supersymmetry (SUSY) with the breaking scale around 1 TeV is a solution to the 
hierarchy problem between the electroweak scale and the Planck/grand unification scale, 
all particles in the standard model have their superpartners with masses not much higher 
than the electroweak scale. These new particles will then be produced at colliders in near 
future, such as the CERN Large Hadron Collider (LHC) and the International Linear 
Collider (ILC). Precision studies of these particles will be possible at these colliders [T|. 

One of the main motivations for these studies is the investigation of the SUSY breaking 
mechanism in the unified theory. The soft SUSY breaking parameters may be extracted 
from precision measurements at future colliders E] > and extrapolated to much higher 
scale jU E] by renormalization group equations. Structure of the soft SUSY breaking at 
higher scale provides an important clue to the SUSY breaking mechanism in the unified 
theory. For example, the unification of three gaugino masses at the same scale as that of 
the gauge couplings is crucial for the SUSY grand unified theory [7j and superstring 
phenomenology jH]. 

For this purpose, we need not only precise measurements of physical parameters such 
as masses and cross-sections, but also precise prediction of the relations between these ob- 
servables and parameters in the lagrangian. Sometimes we have to calculate the relations 
beyond the one-loop order to match expected precision at future experiments. To the 
masses of the particles, for example, two-loop corrections to the top and bottom quarks 

squarks in the first two generations jTUj, and Higgs bosons [TT] have been calculated 
in the framework of the minimal SUSY standard model (MSSM) [12J. 

Here we consider the correction to the mass of the gluino g, the superpartner of the 
gluon, in the MSSM. At the CERN LHC, gluino is, if it is sufficiently light, expected to 
be copiously produced by the strong interaction [T3|. A study [H] has shown that, for the 
SUSY parameter set SPSla given in Ref. [T3] with m~ g ~ 600 GeV, nig may be determined 
to accuracy bm- g = 8 GeV from precision data at the LHC, and even to Srrig = 6.5 GeV 
when combined with the data from the ILC. On the other hand, the one-loop QCD 
correction to the gluino mass [H)J[T7J[TH] is much larger, typically O(10) %, due to large 
coupling constant and large SU(3) representation of the gluino. One naively expect the 
two-loop correction to m g is 0(1) %, similar order to the experimental uncertainty. It 
is therefore important to examine whether higher-order corrections to the gluino mass is 
relevant in extracting SUSY breaking terms at future precision experiments. In addition, 
these corrections might be relevant to the uncertainty ^Hl in the calculation of the SUSY 
particle masses by the computer codes [2H] for given values of the soft SUSY breaking 
terms at the unification scale, at the same order as other two-loop mass corrections and 
the three-loop contributions j2Ij to the running of parameters between the unification 
scale and the scale of the SUSY particle masses. 

In this paper, we calculate the pole mass of the gluino as a function of the lagrangian 
parameters, including SUSY QCD correction to O(a^), by diagrammatic calculation. For 
simplicity, we ignore quark masses and squark mixings in the loops. This approximation 
is valid for the case where the gluino and squarks are sufficiently heavier than the quarks. 
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We also assume degenerate mass rriq for squarks. The effects of the quark masses and 
left-right mixings of squarks to the gluino mass correction will be briefly commented later. 

The pole mass m g of the gluino, which is defined in terms of the complex pole s p = 
(m g — iTg/2) 2 of the gluino propagator, is given at the two-loop order by 

M 3 -E M (s P ) 

= M 3 - Re[M 3 S«(M 3 2 ) + s[f(M 3 2 )] - Re^E^Mf) + 4?(M 3 2 )] 

+Re [{M 3 S«(M 3 2 ) + E$(M 3 2 )}{E«(M 3 2 ) + 2M 3 2 S«(M 3 2 ) + 2M 3 s£ } (M 3 2 )} 

= M 3 + 5mf + 5mf . (1) 

Here M 3 is the running tree-level gluino mass in the lagrangian. E^ M and S^m are the 
one-loop and two-loop parts of the gluino self energy 

iE(p)=i(E x (p 2 )^ + E M (p 2 )), (2) 

respectively. The dot in Eq. denotes the derivative with respect to the external 
momentum squared p 2 . 

The SUSY QCD contribution to E(p) is generated by loops with the gluino, gluon, 
quarks, and squarks. Parameters in the lagrangian are renormalized in the DR scheme 
|22j . which is the DR scheme (dimensional reduction [23] with modified minimal subtrac- 
tion) with additional finite counterterms for squark masses to remove the dependence on 
the mass m e of the e-scalar [21]. 

The one-loop correction 5m~p in our approximation of massless quarks and degenerate 
squarks is [El HE] 



5mf = ^M 3 U - 6 log ^ j + ^-N.TMB^Mi, 0, m 9 ~), (3) 



where Cy — 3, Tp = 1/2, and N q = 6 is the number of quarks. Parameters (a s , M 3 , rriq.) 
in Eq. (JHJ) are the DR running ones at the renormalization scale Q. Bi(p 2 , mi, 1712) is the 
finite parts of the one-loop Passarino-Veltman function [23] , defined by (D = 4 — 2e) 

R( 2 n j Q 2e f d D k k^p 1 

Its explicit form for m\ = is 

(m 2 - - M 2 ) 2 m 2 - - Mo 2 1 m\ m? 

The two-loop 0(a 2 ) correction dm^ consists of two parts, <5m~ 2 ' ) = ^m- 2 ' 1 -* + 5m g 2 ' 2 \ 
(2 1) 

where 5m\ ' is the contribution of the diagrams with only gluons and gluinos, while 
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Figure 1: Two- loop 0(a 2 ) contributions to the gluino self energy, without quark and squark 
propagators. 

<5m~ 2,2 ^ is the remaining contribution with quark and squark loops. Two-loop self energy 
diagrams for these contributions are shown in Fig. Q and Fig. [21 respectively. In these 
figures, the wavy line, solid line without an arrow, solid line with an arrow, and dashed 
line with an arrow represent the gluon, gluino, quark, and squark, respectively. 

The contribution 5m~ 2 ' 1 ' ) is obtained by applying the formula of the 0(a 2 ) QCD cor- 
rection to the quark masses in the DR scheme [2E] to a fermion in SU(3) adjoint repre- 
sentation. The result is 



5m 



(2,1) 



V 4tt 



M 3 ( -48 log -^ + 36 log 



2 -f + 26 + 5tt 2 



4tt 2 log2 + 6C 3 



(6) 



where ( 3 = Y,n=i n — 1-202. We have verified Eq. (|b 
diagrams. At Q — M 3 , the correction © is 5m~ ' /M 3 ~ 



) by explicit calculation of the 
31(a s /n) 2 ~ 0.03. 

(2 2) 

The contribution 5m~ ' including quark and squark loops is calculated as follows: 
The Feynman integrals are decomposed into basic scalar integrals given in Refs. [23 [2E] > 
with the help of the integration by parts technique [23 EH1 CZUJ - The resulting formulas 
are then numerically evaluated by the package TSIL [3T] for two-loop integrals. We have 
analytically checked that the Q dependences of <5m~ 2 ' 1 ' ) in Eq. (JHJ) and 5m g 2 ' 2 ^ are consistent 
with the two- loop renormalization group equation of M 3 . 

(2 2) 

The explicit form of om~ ' is rather long and will be presented elsewhere. Here we 
just show, for reference, the form in the limit of rriq ^> M 3 , obtained by the heavy mass 
expansion technique [531 : 



,(2,2) 

a 2 M- 



5mr* J (m s > M 3j 
72 log 2 



(4tt) 



g+ 2421og^ + log^ 



54 - 288 log -± 



172 + -7T 2 
3 
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Figure 2: Two-loop 0(a 2 s ) contributions to the gluino self energy, with quark and squark 
propagators. Other diagrams obtained by charge conjugation are not shown. 

a*M 3 / 2 M 3 52 M 3 37 4 2 \ 

+ t^ n « CvTf {- 8 log -q + t log "o" - t - r ) ■ (7) 

The last line of Eq. ( |7| ). which is independent of m^, comes from the diagram (a) in 
Fig. |21 We have checked that the rriq dependence of Eq. flj|) is consistent with the two- 
loop running of the gluino mass in the effective theory where squarks are integrated out 

We then present numerical results of the mass correction to the gluino, for the running 
tree-level mass M 3 (M 3 ) = 580 GeV which is close to the values in the SPSla point. The 
strong coupling constant is determined by ot s (mz) = 0.12, which is the running parameter 
within the standard model. 

In Fig. |3J we show the dependence of the one-loop pole mass m~ = M 3 + 8m~ and 

two-loop pole mass rn? = M 3 + 8mp + 5m^ on the renormalization scale Q. We give 

running parameters M 3 (Q ) = 580 GeV and mq(Qo) = 800 GeV at Q = 580 GeV, 

and evolve them and a s to a given Q by 0{a 2 s ) renormalization group equations. For 

reference, the tree-level mass M 3 (Q) decreases from 589 GeV at Q = 400 GeV to 559 

GeV at Q = 1400 GeV. We see that the Q dependence slightly improves as one includes 

(2) (2) 
8m g . One should however note that £m~ is much larger than the Q-dependence of the 

one-loop result mp. It clearly shows that the latter is not an adequate estimate of the 
higher-order contribution to the pole mass. This property has already been observed for 
the 0(a 2 )- correct ion to the squark mass jlOj . 
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Figure 3: Dependence of the one-loop (dashed) and two-loop (solid) pole masses of the gluino 
on the renormalization scale Q. Mass parameters are (M 3 ,m,j) = (580,800) GeV at Q = 580 
GeV. 



In Fig. |4J we compare and as functions of the running squark mass rrig(Q = 
M 3 ), for M 3 (M 3 ) = 580 GeV. Here the renormalization scale is fixed at Q = 580 GeV. The 

(2) 

two-loop correction 5m g is positive and in the range of 8 — 15 GeV for rriq = 400 — 1600 
GeV. This correction is 0(1 — 2) % of the one- loop result rrvp and somewhat larger than 
the expected uncertainty (8.0 GeV to 6.5 GeV) in the gluino mass determination at future 
colliders 14J. For rriq > M 3 , 5m g 2 ^ increases with rriq, as suggested by the rriq 3> M 3 limit 

We also show the result of the modified one-loop formula of the pole mass 



,(2) 



m~ (mod.) 



M 3 (Q) + 



C v a s {Q ) 

47T 



-m- g 



TTlg 



+ ^^NgTrrngBiiml 0, m 5 (pole)), 



7T 



where running masses in Eq. (jSJ) are replaced by the pole masses. Eq. (JHJ) includes higher- 
order corrections by one-loop renormalization group equations ^H] • In Fig. HJ it is seen 
that the modification (JHJ) does not work for the inclusion of the two-loop correction 5rri~^ . 
It would be a useful task to find other modifications of the one-loop formula which incor- 
porate leading part of 5m!?\ 

We finally comment on the effects of the quark masses m q and left-right mixings of 
squarks, both of which are ignored here, to the gluino mass correction. Since these pa- 
rameters break the SU(2)xU(l) gauge symmetry, their contributions to m g should be 
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Figure 4: The pole masses of the gluino at the one-loop (dashed), two-loop (solid), and modified 
one-loop mass Eq. (JHJ) (dot-dashed), for the tree-level mass M 3 (M 3 ) = 580 GeV, as functions of 
mg(M 3 ). 

suppressed by factors m^/m? or m^/m? compared to the SU(2) xU(l)-symmetric contri- 
bution shown in this paper. We therefore expect that, when the gluino and squarks are 
sufficiently heavy, their effects on the 0(a 2 ) correction 5m^ are numerically irrelevant 
for future precision studies of the SUSY particles. However, these effects might become 
relevant for SUSY parameter sets with relatively light gluino or squarks. We will present 
the complete result of the 0(a 2 ) mass correction including these effects, as well as the 
0(a s h 2 ) contributions involving quark-Higgs Yukawa couplings h g , elsewhere. 

In conclusion, we have calculated the two-loop SUSY QCD contribution to the gluino 
pole mass by diagrammatic method, for the gluino and squarks sufficiently heavier than 
the quarks. The 0(a 2 ) correction to the gluino mass has been shown to be typically 
1 — 2 %. For the case of M 3 (M 3 ) = 580 GeV, this correction is similar to, or larger 
than, the expected uncertainty in the mass determination from precision measurements 
at future colliders. This correction would affect the extraction of M 3 from experimental 
data and, since M 3 contributes to the running of many soft SUSY breaking parameters, 
the determination of the SUSY breaking at the unification scale. 

Acknowledgements: The author thanks Mihoko Nojiri for discussion from which this 
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